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Abstract 

This short note aims at (re) proving that the symmetrically normalized graph Laplacian 
L = Id— D~^/'^W D'^/"^ (from a graph defined from a Gaussian weighting kernel on a sampled 
smooth manifold) converges towards the continuous Manifold Laplacian when the sampling 
become infinitely dense. The convergence rate with respect to the number of samples N is 
0(1/7V). 

There exist discrete operators which are the equivalent of the gradient and the divergence 
operators defined on continuous manifold. They share with them some common properties and 
they converge also to their continuous counterparts for a sufficiently fine sampling of the underlying 
manifold. 

The first one relies on the definition of edge derivative [I112]- For a smooth function / : y — M, 
the edge derivative of / on u G y along the edge e = {u,v) £ E formed by the connected vertices 
u and V &V reads 



where d{u) = X^^gy w{u, v) is the degree of the vertex u. From this relation, we have obviously 

Ve/(n) = -Vef{v). 

The gradient of / is then defined globally as the vector field V/ : — )• M defined on the edge 
set E as V/(e) = Vef{u) for e = {u,v) G E. This gradient can be represented as a linear operator 
V e R^x^, so that 

V/(n) = {Vf{u,v) : veV} e R"", 



corresponds to the gradient oi f on u G V seen as a vector of M . The norm of this object on each 
u £ V is defined naturally as 



l|V/(n)f = ^\Vf{u,vr = ^\Vfiu,v)\' 



v<=V 



The scalar products {f,g) = Y2uev f(''^)9('^) i^^G) = YleeE between two real 

functions / and g in the Hilbert space £'^{V) = {h : "^u^y \h{u)\'^ < 00} and two vector fields F 
and G in £'^{E) = {H : Y2eeE < °°}' defined then the adjoint of the gradient, i.e. the graph 
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divergence V*. Indeed, applying the relation (VqjF) = {g, —'V*F) valid in the continuous domain 
[3], we get 

The graph Laplacian of a function / is then defined as 

Using matrix notations, this operator corresponds actually to the common (symmetric) normalized 
Laplacian defined on graph, i.e. 

A = D-^/^WD-^/"^ - Id, 

with D G R^^^ {N = #V) a diagonal matrix such that Duu = d{u) and W G M^^^ the weight 
matrix with Wuv = w{u,v). 

Interestingly, the graph Laplacian converges to the continuous Laplace-Beltrami operator [3] on 
the manifold underlying the graph definition. 

Lemma 1. // the vertex set V = {vi, ■■■ ,vn} C Ai consists of N points taken uniformly and 
independently at random on a m- dimensional compact manifold M embedded in M", then, for 
any smooth function / : F — )• M and for a weighting function iv{u,v) = exp{— ||m — u|p/2e}, the 
normalized graph Laplacian defined by the graph G = (y, E = V xV^w) satisfies 



1 hm M{n) = \AMf{u) + 0{e 



1/2N 



with A_\4 the Laplace-Beltrami operator defined on M . 

Proof. We follow a similar development to the one given in [5] . With the hypothesis of the Lemma, 
for any function h G I'^iV), we have 

\Wh]{u) = '^w{u,Vj)h{vj) = '^eyi^{-\\u - VjW'^ /2e}h{vj) + h{u). 

Since Yj = w{u, vj) h{vj) are iid for j ^ i, by the law of large numbers we have 
'^w{u,Vj) h{vj) ~ (iV - l)E^[exp{n - 



_ N-1 
volM 



/ exp{-||u-y|p/2e}/i(?/) dy^iy, 
JM 



where the integral is performed on the manifold with the local infinitesimal volume element dxy. 
The relation (2.9) in [5] (or Eq. (10) in [6]) explains that 



(2^)^/^ j^exp{-\\u-yf/2e}h{y) ^mV = h{u) + ^[E{u)h{u) + AmKu)] + 0{e 
where E{u) = \S{u) and S is the scalar curvature of jl]. Therefore, 



m]{u) = [h{u) + ^,[E{u)h{u) + AMh{u)] + 0(6^/2)] + h{u) 

= ^""'S^""' Hu) + m^)h{u) + AMKu)]+0{e'/\l/{m"^/'))], (1) 
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where the notation 0{^,i') means |0(/i, < C/i + Dv for two positive values <^ 1. Taking 
/i = 1, we get then 

d{v,) = ^^^^^M^ [1 + f E{v,) + 0{e'l\ l/{Ne-/^))] . (2) 

Therefore, 



Taking now = f{t)g{t) ioi t ^ M with 5r(t) = —. — == in dH), we get 

= [g{u)f{u){l + ^E{u)) + ^AmHu) + 0(e3/2, l/(iVe-/2))] , 



(Ar-l)(27re)'"/2 



However, 

AmHu) = f{u)AMg{u) + 2 (VA4/(n), V>i5(u))T„>i + g{u)AMf{u), 

where the scalar product (•,-)t„A4 occurs in the tangent plane T^M ~ M'^ of on n with the 
metric of the manifold |3]. For any function s on A^, the gradient V m^{'^) ^ is composed of 

the derivative of s according to a local system of coordinates (or chart) isomorphic to M'^. 

From the definition of it is clear that any differential operator Dj^^^ of g with respect to this 
local system in TuM. is of order Dj^g{u) = 0(e). Consequently, 

^AmHu) = ^g{u)AMf{u) + 0(e2), 

that provides the final result, 

Af{u)= g\u)f{u){l+^^E{u)) + fA^,/(n) - f{u) + 0{^l\l/ {Ne"''^)) 
= ^AMfiu) + 0(e3/2,l/(Are™/2)). 

□ 
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